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Abstract 

Ohmic currents induced prior to decoupling are investigated in a standard transport 
model accounting both for the expansion of the background geometry as well as of its rel- 
ativistic inhomogeneities. The relative balance of the Ohmic electric fields in comparison 
with the Hall and thermoelectric contributions is specifically addressed. The impact of the 
Ohmic currents on the evolution of curvature perturbations is discussed numerically and it 
is shown to depend explicitly upon the evolution of the conductivity. 



Prior to photon decoupling the plasma is electrically neutral, the static (Coulomb) po- 
tential is exponentially suppressed beyond the Debye length (see, e.g. ^) while the concen- 
tration of the electric charges is 10~^° times smaller than the concentration of the photons 
(see e.g. [2]). This effect would naively seem to increase the role of the Hall and thermo- 
electric terms whose magnitude is inversely proportional to the charge concentration [5]. 
Ohmic electric fields might also be induced because of the presence of large-scale magnetic 
fields. The value of the conductivity is then crucial for determining the magnetic and electric 
diffusivity scales. The aim of the present paper is to clarify the situation and investigate 
more quantitatively the different contributions responsible of Ohmic current^] especially in 
the light of the ongoing attempt of a consistent inclusion of large-scale magnetic fields in 
the calculation of the Cosmic Microwave Background (CMB) observables HE]- Consider, 
to begin with, the Vlasov-Landau system of equations for electrons and ionqj 

^ + tT- V.^/e,i Te[E + v^B]- V,-/e,i = (^)^^^^, (1) 

where E = o?E and B = a^B are, respectively, the comoving electric and magnetic fields; 
V = qj ^f^^^rm^ is the velocity and q is the comoving three-momentum. In the ultra- 
relativistic limit (i.e. q ^ ma), v = q/\q\ and, therefore, Eq. ([T]) is invariant under a 
Weyl rescaling of the geometry g^j^^: this boils down to the conclusion that, absent the 
relativistic fluctuations of the geometry (which will be introduced in a moment) the Vlasov- 
Landau system has the same form it would have in fiat space-time provided the underlying 
background geometry is spatially fiat. Conversely, q = mav when the given species are non- 
relativistic; Weyl invariance is then broken by the masses of the electrons and of the ions 
(i.e., respectively, m^. and m\ ~ mp). In terms of the distribution functions of Eq. ([1]), the 
evolution equations of the electromagnetic fields are given by: 

V ■ E = 47re 1 d\[f,{x, v, r) - /e(f , v,t)], V ■ 5 = 0, (2) 

V X E + 5' = 0, V X B - E' = Ane J d\v [fi{x, v, r) - /e(f , v, r)], (3) 

where the prime denotes a derivation with respect to the conformal time coordinate r. The 
evolution equations of the comoving concentrations of electrons and ions (i.e. respectively 
He and rii) can be written, in explicit terms, as^ 

On. OtI' -» 

^ + V ■ (netJe) - 3ne^' = 0, ^ + ^ ■ (^i^i) - = 0, (4) 



^For cold plasmas in the laboratory the so-called Ohm law (see, e.g. [T]) is a relation (often non-linear) 
involving the total current the electric and magnetic fields, the electron and ion pressures, the bulk velocity 
of the plasma. In curved backgrounds, on top of the previous quantities, one needs to take into account the 
effects due to the expansion and to the large-scale fluctuations of the geometry. The resulting Ohm law is 
then, effectively, an evolution equation for the Ohmic current. 

^The conformal time coordinate will be denoted by r and the geometry will be assumed to be conformally 
flat, i.e. gfj^i, = a'^{T)rj^^, where r/^^ = diag(l, —1, —1, —1) is the Minkowski metric. 

■^Throughout the paper the physical quantities will be denoted by a tilde while the comoving quantities 
will appear without the tilde. For instance the (comoving) concentrations, energy densities and pressures of 
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where if) is the (scalar) fluctuation of spatial components of the metric in the longitudinal 
gauge [To] defined by the conditions 5s5'oo = 2a^0, dsQij = 2a?ip5ij. Introducing the global 
charge and the total current, i.e. 

Pq = e(ni-ne), J = e{niVi - rieVe) , (5) 

the difference of the two equations reported in Eq. (jl]) implies that + V ■ J — 3ip'pq = 0. 
Using Eq. the relevant Maxwell equations become V ■ E = 47rpq and V x 5 = AnJ + E'. 
Recalling that the pre-decoupling plasma is globally neutral, i.e. ric = rii = rji^n^ where is 
the comoving concentration of photons and rj^ is the ratio between the baryonic concentration 
and the photon concentration, i.e. r]^ = 6.29 x 10"^° (/igfibo/0.02273)[T^o/(2.725 K)]-^ where 
Obo is the critical fraction of baryons and T^o is the CMB temperature. The fiducial values 
of the cosmological parameters employed to illustrate the present estimates correspond to 
the best fit of the WMAP Syr data alone [6l[7j. The conductivity (and the related mobility) 
can be computed in the customary framework of the Krook model [H [9] which holds for 
weakly ionized plasmas, and, with some numerical differences, also in the fully ionized case. 
The collision terms of Eq. ([T]) can then be written as 



. ; rei(/e-/e), hf =rie(/i-/i), (6) 

or / coll \or / coll 

where Fei and Fie are the collision rates of electrons and ions and where /g ; are two 
Maxwelhan distributions, i.e. f{v) = [ma/(27rT)]^/^ exp [— mau^ a/(2T)]. The induced elec- 
tric field slightly perturb the Maxwellian distributions and, therefore, the explicit form of 
the conductivity can be derived from Eq. ([3]) by following exactly the same steps of the 
standard calculation plasma calculation (see e.g. [H [3]) with the important difference that, 
because of the breaking of Weyl invariance, the scale factors a(r) appear ubiquitously: 



CT = — ^^W^[lnAc(T)]-^ = 4.35 x 10"^ eV ^^^^ 
SvTA/Se^V m,a^ V2.725Ky V0.1326y 



871 V3> 

-1/2 



9 TIT.,. T^\3/2 ./^2^y/2 ^ 

V2.7r"" 

Ac(T) = Af—V^' = 1-102 xlO^f^^V , (7) 
^ 2e^\7moJ V0.02273/ ' ^ ' 

where Hq is the common value of the (comoving) electron and ion concentrations; Ac is the 
argument of the Coulomb logarithm and ^mo is the critical fraction of matter in the ACDM 
model. We are now interested in the evolution equation of the Ohmic current J whose 
explicit form can be derived by combining the governing equations for electrons and ions: 

+ nvi = —[E + Vi X B] - V<P - ^ + ar,A{ve - v,) + -^aT.Jv^ - v,), (8) 
Pitt api Pi 3 Pi 

Ve ' + nve = -[E + tJe X fi] - V0 — + aT,i{vi - V,) + -^aT,^{v^ - V,). (9) 

Pea ape 3 pe 

electrons and ions will be denoted, respectively, by rtc,i — a'^'^c.i, Pc,i — o? Pc\ and by po, i = ?^c, iTc, i = a^pcA- 
Similarly, for photons, — (7r^/15)r^ = a^p-f. When needed these two notations will be employed without 
further explanations. 
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By taking the difference of Eq. ([S]) (multiplied by erij) and of Eq. 
the following equation can be obtained: 



(multiplied by erie 



dJ 
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where the plasma frequencies and the baryonic velocity have been introduced: 
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me + mj. 



(10) 



The evolution of Wb is coupled to the velocity of the photons and it is obtained by summing 
up (instead of subtracting) Eq. ^ (multiplied by mi) and Eq. (I9j) (multiplied by me): 



Jx B 
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aVb(l + rric/mi) 

V4> + aV^e{v\, — Vj 



3pb 



(12) 
(13) 



Eq. ffTOj) can be expanded in power in powers of (me/mi). Recall that pb = mehe + rriihi and 
that, by global neutrality, ni = ne = no where no = '/7b^7 where r^b is the ratio between the 
baryonic concentration and the photon concentration already introduced after Eq. ([5]). The 
result of this double expansion implies, from Eq. (fTOj) . 
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The terms J' and TiJ are comparable in magnitude and are both smaller than Fie and Fe^, 
i.e. l-LJ ~ J' < (4/3)(p^/me)Fe^ < aFie. While it is important to solve the evolution 
of J during all the pre- decoupling regime, the previous chain of inequalities implies that, 
asymptotically, the form of Eq. f|T^ is dominated by the term containing Fie. At the right- 
hand side the term containing [y^ — v^) can be estimated by subtracting Eqs. ( fT2l) and ( fT3l) . 
The difference (wb ~ ^7) is driven exponentially to zero at a rate controlled by aF^e(l + -^b ^) 
where i?b = (3/4)pb/p^- The asymptotic form of the Ohm's law can be written, for large 
conformal times as 

,2 
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UJ 



pc 



47r[aFie + (4/3)(p^/me)FeT,] 



E + VhX B + 



eno 



J X B 



(15) 
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The term containing the gradient of the electron pressure is the curved-space counterpart of 
the thermoelectric term [3] while the term proportional to the vector product of the current 
and of the magnetic field is the curved-space counterpart of the Hall term. The displacement 
current can be neglected in comparison with the Ohmic current (i.e. 4nJ E') provided 
the left hand side of Eq. (fn|) is subleading in comparison with the induced electric field 
(i.e. J' <^ ujp^^E). The latter requirement demands, after derivation with respect to the 
conformal time r, the fulfillment of the condition J" <^ ^pc^'- one-fluid description 
correctly captures the dynamics in the low-frequency branch of the spectrum of plasma 
excitations, i.e. uj -C cUpe- If the thermoelectric and Hall terms are neglected, then the 
electromagnetic fields obey the following pair of equations, i.e. 

OB ^ , -X 1 9 - dE d , 1 , - 

— = V X 4 X i? + -—V^B, — = -— 4 xB) + -—V'E, 16 

or Ana or or Ana 

where a = uj"^^/ (47rrei) is given by Eq. ([7]); Eq. f|T6|) implies that wavenumbers k"^ > = aH 
are dissipated because of the finite value of the conductivity. The explicit value of the 
diffusivity scale k^- is given by 

^ = 2,177 X 10-- f^V'7^i5M^)="V»/*. a = (17) 
T V10.75/ V0.1326/ ' Oeq ^ ^ 



In comoving temperature units, the Hubble wavenumber is 7i/T = 1.047 x lO^^^y (7p/10.75 
(where Qp is the effective number of relativistic degrees of freedom); as expected not only 
7i <^k^ but also Ti -C fcoebye — 5.26 x 10~^T where /cDebye is the wavenumber corresponding 
to ADebye = \Jt / {SnuQe^) i.e. the screening length of the Coulomb potential between two 
charges in the plasma. The diffusivity scale, on the contrary, sets a (lower) limit in the 
coherence scale of the Ohmic fields. The dominance of the drift term (i.e. I^b x -B|) over the 
thermoelectric and Hall terms demands, from Eq. f[T^ . the fulfillment of the following pair 
of relations: 



In Fig. [T] (plot at the left) the contribution of the Hall and thermoelectric terms to the 
Ohmic current is illustrated. For comparison the diffusivity, Debye and Hubble scales are 
also reported. The shaded area denotes the region where the conditions of Eq. f|T8l) are 
approximately fulfilled, i.e. for typical amplitudes of the comoving magnetic field in the 
range 10~^nG < B < lO^nG. Because of the value of the charge concentration, the Hall 
contribution dominates over the thermoelectric term of the electrons 

\JxB\ B^ 9 / 5 \ V Bn^o \ , ^ 

' ' - 8.19 X 10^ — L ° • (19) 



|Vpc| 47™oT VnG/ V 0.02273/ 

When the plasma contains a magnetic field whose Fourier modes are stochastically dis- 
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Figure 1: Relative contribution of Hall, thermoelectric and drift terms in the Ohmic current 
(plot at the left). Contribution of Ohmic currents to the evolution of curvature perturbations 
(plot at the right). In the left plot (as well as in Fig. [2]) ti sets the equality time scale (more 
precisely Teq = (v^- l)ri). 




Figure 2: The contribution of the Ohmic fields to the curvature perturbations (plot at the 
left) and to the metric fluctuations (plot at the right). 



tributed with power spectrum Pb(^), the asymptotic form of the Ohm law ( |T5i) . within the 
shaded region of the parameter space of Fig. [U induces an effective electric field 
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(20) 



where Pij{k) = {6ij — kikj/\k\^) is the transverse projector. The stochastic electromagnetic 
fields as well as the induced Ohmic currents, being inhomogeneous, affect the curvature 
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perturbations whose evolution, on the absence of non-adiabatic pressure fluctuations, is 
given by [H [5] : 

E-J 7^(5sPb(34-1) , [?,clgW + gl{^cl-2)] {V ■ v,) 

^ 3a4(pt+pt) 3(pt+pt) 3(pt+pt) gla^ + gl 3 ' ^ ^ 

where 5sPb = -B^/(87ra^) and 5sPe = -E'^/(87ra^); moreover 

Note that (pt + pt)vt = Y^aips. + Pa)Va. is the total velocity field of the plasma including the 
contribution of cold dark matter particles, neutrinos, electrons, ions and photons. When the 
Universe contains matter, radiation and dark energy the total barotropic index Wt and the 
total sound speed can be written, respectively, as 

* Pt 3(ai + a4 + aai)' pi * 3(1 + u-t) 9a ' ^ ^ 

where, as already stressed, a = aja^^] furthermore Pt = Pr + Pm + Pa and 

Q,^ _00/l«C1( U "^'J \ I U A ^ ^, _ Q1 OK 1 O I U MU \ 

'•''eq 



::^ = 2246.81 ^3)::^ = 3195.I8 pi::^ . (24) 

Oen V0.1326/ V0.3835/ ' V0.1326/ ^ ' 



The evolution of C, as well as the evolution of if) can be integrated and the results are reported 
in Fig. [T] (plot at the left) and in Fig. [21 Both in Fig. [T] and [2] f^B = SsPB/iSnp^) and 

= 4pE/(8vrp^). In Fig. [T] (plot at the left) the pure magnetic contribution is compared 
to the total Ohmic contribution in units Qe = 1- In Fig- [2] the contributions of ■ J is 
compared to the other terms arising in the evolution of ( and tp. In the baryon rest frame 
the Ohmic contribution is suppressed as (/c^/a^) ~ (9(10^^°) for length-scales larger than 
the Hubble radius (notice that, indeed, in the left plot of Fig. [2] has been rescaled by 
a factor 10^° to make the two contribution visually comparable on a linear scale). The 
latter result is compared with the suppression experience by ■ J which is of the order of 
(k'^/kl) ~ O(10~^^). In Fig. [2] (plot at the right) the contribution of E ■ J (rescaled by a 
factor 10^^) is compared with the magnetic contribution as it arises in Eq. (12T]) . The Ohmic 
contribution is dominated by the drift term which vanishes in the baryon rest frame and 
which is subleading over typical length-scales larger than the Hubble radius. An interesting 
byproduct of this study is the derivation of a consistent evolution equation for the Ohmic 
current. The latter result improves on the usual approximations posited in the Boltzmann 
integrators accounting for the effects of large-scale magnetic fields [H |3] on CMB observables. 

NQL wishes to thank the CERN physics department and, in particular. Prof. L. Alvarez- 
Gaume for kind hospitality and financial support. 



''We shall denote by C the density contrast on uniform curvature hypersurfaces (see, e. g. [TT]) while TZ 
is the curvature perturbation on comoving orthogonal hypersurfaces. The two quantities are connected by 



the Hamiltonian constraint, i.e. Eq. ((22|) . 
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